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ABSTRACT 

We develop a method for predicting the yield of transiting planets from a photometric survey given the pa- 
rameters of the survey (nights observed, bandpass, exposure time, telescope aperture, locations of the target 
fields, observational conditions, and detector characteristics), as well as the underlying planet properties (fre- 
quency, period and radius distributions). Using our updated understanding of transit surveys provided by the 
experiences of the survey teams, we account for those factors that have proven to have the greatest effect on the 
survey yields. Specifically, we include the effects of the surveys' window functions, adopt revised estimates of 
the giant planet frequency, account for the number and distribution of main-sequence stars in the survey fields, 
and include the effects of Galactic structure and interstellar extinction. We approximate the detectability of a 
planetary transit using a signal-to-noise ratio (S/N) formulation. We argue that our choice of detection criterion 
is the most uncertain input to our predictions, and has the largest effect on the resulting planet yield. Thus 
drawing robust inferences about the frequency of planets from transit surveys will require that the survey teams 
impose and report objective, systematic, and quantifiable detection criteria. Nevertheless, with reasonable 
choices for the minimum S/N, we calculate yields that are generally lower, more accurate, and more realistic 
than previous predictions. As examples, we apply our method to the Trans-Atlantic Exoplanet Survey, the 
XO survey, and the Kepler mission. We discuss red noise and its possible effects on planetary detections. We 
conclude with estimates of the expected detection rates for future wide-angle synoptic surveys. 

Subject headings: methods: numerical — planetary systems — surveys — techniques: photometric 



1. INTRODUCTION 

There are four ways by which extrasolar planets have been 
detected. The first method to unambiguously detect an ex- 
trasolar planet was pulsar timing, which relies on detecting 
periodic variations in the timing of the received radio signal 
that occur as the pulsar orbits about the system's barycen- 
ter. The first system of three planets was fou nd around PSR 
B1257H-12 in 1992 ( Wolszczan & Frai]|ll992j), followed by a 
single planet around PSR B1620-26 (iBacker et alJll993h . Al- 
though rare, the pulsar planets are some of the lowest mass 
extrasolar planets known: PSR B1257H-12a is about twice the 
mass of the Moon. 

The second way to find extrasolar planets is through radial 
velocities (RV), which uses the Doppler shift of observed stel- 
lar spectra to look for periodic variations in the target star's ra- 
dial velocity. By estimating the mass of primary star, the ob- 
served radial velocity curve and velocity semi-amplitude can 
then be used to directly calculate the inclination-dependent 
mass (MpSin;) of the companion object. To date, RV sur- 
veys have detected more than 230 planets around other stars, 
making it the most successful method of extrasolar planet de- 
tection. While the large number of detected systems having 
an unseen companion with a mass on the order of lMy„p sin / 
statistically ensures that the majority of these are planetary 
bodies, the RV surveys are unable to provide more informa- 
tion than the minimum masses, periods, eccentricities, and the 
semi-major axes of the planets. 

RV surveys also have a limited ability to detect planets 
much smaller than a few Earth masses. The state of the art 
in RV surveys is the High Accuracy Radial Velocity Planet 
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Searcher (HARPS) spectrometer at the La Silla Observatory 
in Chile, which is capable of radial velocity measurements 
with precisions better than 1 ms"' for extended periods of 
time (ILoviset all 120061) . HARPS is therefore able to detect 
planets with masses on the order of 3 to A-M(^ in relatively 
short period orbits. Unfortunately, planets closer to an Earth 
mass will be increasingly difficult to detect since intrinsic stel- 
lar variability, in the form of acoustic oscillation modes and 
granulations on the photosphere, makes more precise spectro- 
scopic radial velocity measurements harder to acquire. How- 
ever, it may be possible to surmount this obstacle (as in the 
case of iLovis et al.ll2006l) through the selection of stars with 
"quiet" photospheres and long integration times which serve 
to average out the stellar variability. 

Gravitational microlensing is another technique for detect- 
ing extrasolar planets. Microlensing of a star occurs when 
a a star passes near the line of sight of the observer to an- 
other background star. The gravity of the foreground star 
acts as a lens on the light emitted by the background star, 
which causes the star in the background to become momen- 
tarily brighter as more light is directed towards the observer. 
Planetary companions to the lens star can further magnify 
the background star, and create short-term perturbations to 
the microlensing light curve. To date, six planets have been 
detec t ed using microlensing (iBond etal. l2004t lUdalski et ak 
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| 2005[ iBeaulieu et all l2006riGould et al.ll2006at iGaudi et al 
■2008"). Unfortunately, the one-shot nature of microlensing ob- 
servations means that information about systems discovered 
this way is generally sparser than that available for RV sys- 
tems. Therefore, microlensing is most useful in determining 
the general statistical properties of extrasolar planets, such as 
their frequency and distribution, and not the detailed proper- 
ties of the planetary systems. 

Planetary transits are a fourth method by which extrasolar 
planets have been discovered, and the one that provides the 
most complete set of information about the planetary system. 
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Only planets with very specific orbital characteristics have a 
transit visible from Earth, because the orbital plane has to be 
aligned to within a few degrees of the line of sight. Therefore 
transiting planets are rare. Nevertheless, a transiting extraso- 
lar planet offers the opportunity to determine the mass of the 
planet (when combined with RV measurements) since the in- 
clination is now measurable, as well as the planetary radius, 
the density, the composition of the planetary atmosphere, the 
ther mal emission from the planet, and many other properties 
(see ICharbonneau et alj (12007) for a review). Additionally, 
and unlike RV surveys, transiting planets should be readily 
detectable down to IR(q and beyond, even for relatively long 
periods. 

Having accurate predictions of the number of detectable 
transiting planets is immediately important for the evaluation 
and design of current and future transit surveys. For the cur- 
rent surveys, predictions allow the operators to judge how 
efficient are their data-reduction and transit detection algo- 
rithms. Future surveys can use the general prediction method 
that we describe here to optimize their observing set-ups and 
strategies. More generally, such predictions allow us to test 
different statistical models of extrasolar planet distributions. 
Specifically, as more transiting planets are discovered and 
characterized, predictions relying on incorrect statistics of ex- 
trasolar planet properties will increasingly diverge from the 
observed set of distributions. 

Using straightforward estimates it appears that observing a 
planetary transit should not be too difficult, presuming that 
one observes a sufficient number of stars with the requi- 
site precision during a given photometric survey. Specifi- 
cally, if we assume that the probability of a short-period giant 
planet (as an example) transiting the disk of its parent star 
is 10%, and take the results of RV surveys which indicate 
the fr equency of such planets is about 1% (ICumming et alJ 
|2008|) . together with the assumption that typical transit depths 
are also about 1%, the number of detections should be w 
10~-'A^<i%, where N<\% is the number of surveyed stars with 
a photometric precision better than 1 %. 

Unfortunately, this simple and appealing calculation fails. 
Using this estimate, we would expect that the TrES survey, 
which has examined approximately 35,000 stars with better 
than 1 % precision, to have discovered 35 transiting short pe- 
riod planets. But, at the date of this writing, they have found 
four. Indeed, overall only 51 transiting planets have been 
found at this time by photometric surveys specifically de- 
signed to find planets around bright stars'*. This is almost 
one hundred times less than what was originally predicted by 
somewhat more sophisticated estimates (.Horne.2003h . 

Clearly then, there is something amiss with this method 
of estimating transiting planet detections. Several other 
authors have developed more comple x models to predic t 
the expected yields of transit surveys. iPepper et al.l ( l2003h 
examined the potential of all-sky surveys, which was ex- 
panded upon and generalized for photometric searches in 
clusters 'Pepper & Gaudil (l2005h . iGould et all (12006^ ) and 
iR^essin et al. (2007) tested whether the OGLE planet detec- 
tions are statistically consistent with radial velocity planet dis- 
tributions. Brown (2003) was the first to make published es- 
timates of the rate of false positives in transit surveys, and 
iGillon et al] ( l2005h model transit detections to estimate and 
compare the potential of several ground- and space-based sur- 

As of June 2008. See the Extrasolar Planets Encyclopedia at 
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veys. 

As has been recognized by these and other authors, there 
are four primary reasons why the simple way outlined above 
of estimating surveys yields fails. 

First, the frequency of planets in close orbits about their 
parent stars (the planets most likely to show transits) is likely 
lower than RV surveys would indicate. Rece nt examina- 
tions of the results from the OG LE-III field by ' Gould et al] 
(2006b) and iFressin et al.l (l2007h indicate that the frequency 
of short-period Jovian-worlds is on the order of 0.45%, not 
1.2% as is often ass i imed by extr apolating from RV surveys 
(iMarcy et al.ll2005al) . iGould et all (|2006b) point out that most 
spectroscopic planet searches are usually magnitude limited, 
which biases the surveys toward more metal-rich stars, which 
are brighter at fixed color These high metallicity stars are 
expected to have more planets than solar -metallicity stars 
dSantos et al.ll2004tlFischer & Valentill2005h . 

Second, a substantial fraction of the stars within a survey 
field that show better than 1% photometric precision are ei- 
ther giants or early main-sequence stars that are too large 
to ena ble detect able trans it dips from a Jupiter-sized planet 
(Goul d & Morga n 2003; Brownll2003h . 

Third, robust transit detections usually require more than 
one transit in the data. This fact, coupled with the small frac- 
tion of the orbit a planet actually spends in transit, and the 
typical observing losses at single-site locations due to factors 
such as weather, create low window probabilities for the typ- 
ical tran sit survey in the majority of orbital period ranges of 
interest (Ivon Braun & Ciardill2007h . 

Lastly, requiring better than 1% photometric precision in 
the data is not a sufficient condition for the successful detec- 
tion of transits: identifiable transits need to surpass some kind 
of a detection threshold, such as a signal-to-noise ratio (S/N) 
threshold. The S/N of the transit signal depends on several 
factors in addition to the photometric precision of the data, 
such as the depth of the transit and the number of data points 
taken during the transit event. Additionally, ground-based 
photometry typically exhibits substantial auto-correlation in 
the time series data points, on the timescales of the transits 
themselves. This additional red noise, which can come from 
a number of environmental and instrumental sources, sub- 
stantially reduces the statistical power of the data dPont et al] 
l2006h . 

In this paper we describe a method to statistically calculate 
the number of planets that a given transiting planet survey 
should be able to detect. We have done so in a way that maxi- 
mizes the flexibility of our approach, so that it may be applied 
in as many circumstances as possible. Therefore, we have al- 
lowed for survey-specific quantities such as the observation 
bandpass to be defined arbitrarily, and have also allowed for 
our astrophysical assumptions to be easily altered. We ac- 
count for factors such as the low frequency of Jovian plan- 
ets, the increased difficulty of detecting transits at high im- 
pact parameters, as well as variations in stellar density due to 
Galactic structure and extinction due to interstellar dust. We 
also use the present-day mass function to realistically account 
for the number of main-sequence stars that will be in a given 
field of view. We do not treat, or simulate, the distribution 
and effects of giant stars. We approximate the detectability of 
a planetary transit using a S/N formulation. 

This differs from the previously published work in the field 
by virtue of the generalized nature of our approach . Othe r 
authors, such as jPepper et alj (l2003l) iGillon et all (l2005h . 
IGould et alJ (l2006bl) . and lFressin et all (l2007h have restricted 
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themselves to modeling specific transit surveys or types of 
transit surveys. Our approach, which allows for variations 
such as arbitrary telescope parameters, observing bandpasses, 
cadences, and fields of view, is not restricted to any one 
type of transit survey. Indeed, we offer examples of how 
our method may be applied to the four most frequently used 
modes for transit searches. Furthermore, we have conducted 
an extensive literature review of basic astrophysical relations, 
such as the Present Day Mass-Function (PDMF), Galactic 
structure, and the main sequence mass-luminosity relation, 
and included these in our predictions. We do not use more 
complicated (but presum ably more accurate) G alactic models 
like the Besancon model ( Robin & Crezell986 ') because these 
models are less flexible and considerably more time consum- 
ing to utilize. 

We first describe the mathematical formalism with which 
we have chosen to address this problem, and then move on 
to discuss our specific assumptions. As examples, we then 
offer our predictions for several different types of photomet- 
ric surveys: point-and-stare, drift-scanning, space-based, and 
all-sky. For each, we choose a representative survey for that 
particular mode, and compare our predictions to actual survey 
results, or to predictions elsewhere in the literature. We also 
examine the possible effects of red noise on our predictions, 
and conclude by discussing what our predictions imply for the 
general state of photometric transiting planet surveys. 

2. GENERAL FORMALISM 

In the most general sense, we can describe the average num- 
ber of planets that a transit survey should detect as the proba- 
bility of detecting a transit multiplied by the local stellar mass 
function, integrated over mass, distance, and the size of the 
observed field (described by the Galactic coordinates il,b): 



d'Nde, 



-=p^(r,l,b) r cosb 



dn df(Rp,p) 



dRp dp dM dr dl db dM dRp dp 

xPde,{M,r,Rp,p), (!) 

where Pdc,(M,r,Rp,p) is the probability that a given star of 
mass M and distance r orbited by a planet with radius Rp and 
period p will present a detectable transit to the observing set- 
up, j^ "f is the probabihty that a stai" will possess a planet of 

radius Rp and period p. dn/dM is the present day mass func- 
tion in the local solar neighborhood, and p* is the local stellar 
density for the three-dimensional position defined by {r,l,b). 
We use r^ cos b instead of the usual volume element for spher- 
ical coordinates, r^ sin (j>, because b is defined opposite to 0: 
b = 90° occurs at the pole. 

First, we will integrate over both planetary radius and or- 
bital period 
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" df(Rp,p) 
'R„,„i„ Jp^,n dRp dp 
X Pde,(M,r,Rp,p)dRpdp. 



(2) 



To do this, we must first specify the detection probability 
Pdei, as well as the distribution of planet radius and periods 
df{Rp,p)/dRpdp, which we treat in Section 2.2. 

2.1. The Detection Probability Pdet 



The probability that a planet orbiting a given star shows a 
detectable transit may be broken up into two separate prob- 
abilities: the chance that a planet with radius Rp and period 
p around a star of mass M at a distance r will transit with 
a sufficient signal-to-noise ratio to be detected, and the win- 
dow probability that a transit will be visible for the particular 
observing set-up: 

Pde,{M, r, Rp,p) = Ps/NiM, r, Rp , p) P.Up) 0) 

2.1.1. Ps/N 

To calculate the probability that a planetary system will 
show detectable transits, we first select the appropriate sta- 
tistical test to determine whether or not our notional data 
shows a discernible transit. In the world of transiting planet 
searches, many different methods are used either individu- 
ally or together to test the statistical significance of a possi- 
ble transit and alert the researcher to its existence. For ex- 
ample, the Hungarian-made Automated Telescope Network 
(HATNet) uses a combination of the Dip Significance Pa- 
rameter (DSP) and a Box-Least Squares (BLS) power spec- 
trum (amongst others) to test its photometric data. In general 
though, whatever test one uses, it is generally the S/N of a 
transit shape in the observed data is what determines whether 
or not that transit will be identified as an event. 

As the square of the signal-to-noise can be approximated 
by the x^ test, we therefore chose to use the x^ value to de- 
termine if a given transit will be detectable for a set of given 
observation parameters. For a specific transit, the x^ statistic, 
assuming uncorrected white noise, ^ is given by 

2 



X 



■■N,r 



(4) 



Where A^„ is the number of data points observed in a transit 
of fractional depth S and with fractional standard deviation a, 
assuming a boxcar transit shape with no ingress or egress. To 
be detected, the transit must have a x^ higher than a certain 



mmimum x„ 



A — Xmi 



(5) 

The fractional depth of the transit is given by the square of the 
planet-to-star radius ratio 

2 



s = 



(6) 



While this formulation ignores the contributio ns of stel- 
lar limb-darkening to the depth of the transit, Gould et alJ 
(2006b) demonstrate that, to a first order approximation, tran- 
sit detection rates are not affected by limb-darkening. Briefly, 
they point out that the effects of limb-darkening on the transit 
lightcurve will be offset by the effects of the planet's ingress 
and egress from the stellar disk. Therefore, we have not in- 
cluded stellar limb-darkening terms in our calculations. 

To determine a, we assume Poisson statistics, such that for 
a given exposure with an observed number of A^s source pho- 
tons from the target star and A^^ background photons. 



Ns+Nb 



N; 



+ r ■ 

2 scint ' 



(7) 



where crj„„, accounts for scintillation as per lYoungl (Il967h . 
and Ns is defined as 

Ns = exFs.xte.pA. (8) 

' We consider red noise later, in Section 4.5. 
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Here e\ is the efficiency of the overall observing set-up, and 
can take values from to 1. It accounts for photon losses 
in places such as the filter, the mirrors of the telescope, the 
CCDs, and passage through the atmosphere. Fs.\ is the source 
photon flux in the observed bandpass in units of 7 m~^ s~^ . texp 
is the exposure time of each observation, and A is simply the 
light collecting area of the telescope. 

Fs,x, the observed photon flux from the source in a certain 
bandpass, can be defined as a function of the photon luminos- 
ity of the source, the distance to the source, and the effects of 
extinction due to interstellar dust: 



Fs.\ = 






10 



-Ax/2.5 



(9) 



where 10"'*^/^^ accounts for extinction, and will be treated 
later It has the effect of decreasing the number of observed 
photons from more distant stars, particularly in the bluer 
bands. 

^\(M), the photon luminosity of the source star, we set as a 
function of the stellar mass M and the observational bandpass 
A. We derive this quantity in Appendix A, and only note here 
that we approximate the source stars as blackbodies. While 
this is a simplification, we show in Appendix A that it is suf- 
ficient for our purposes. Using ^x{M) in equation ^ yields 
the observed photon flux at a distance r from that same star. 

We thus now have an expression for Ns as a function of the 
stellar mass, distance, and bandpass choice, that we can use 
in our calculation of a. We now only need A^^, the number of 
background photons observed. Similar to equation (O, A^^ is 
defined as 



Nb = Sskv.X^texpA, 



(10) 



where Ssky.x is the photon surface brightness of the sky at the 
particular wavelength A (or, in our case, for a specific band), 
and n is the effective area of the seeing disk. Assuming that 
the point-spread function (PSF) is a Gaussian with a full- 
width half-max (FWHM) of Qfwhm arcseconds, J7 is defined 
by 



'^ - ]j.A^FWHM- 



(11) 



We can now place both of the expressions for Ns and A^^ back 
into equation (Q to arrive at a determination of a, the photo- 
metric uncertainty in the data points. 

In certain cases, Ns and A^^ will be so great that the CCD 
detectors of the telescope will be saturated by the number of 
photons received from a star during the exposure time. In this 
case, further measurements become useless. Mathematically, 
we may describe saturation as occurring when the number of 
photons collected by a single CCD pixel during an exposure 
exceeds the Full-Well Depth (FWD) of that pixel: 



N pixel > NfWD- 

N pixel is related to Ns and Nb as 



(12) 
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(13) 



where 8p,v is the angular size of an individual pixel. 

The only part of equation (|4]i that remains to be defined 
is the number of data points observed while the system is in 
transit, N,,-- Assuming a circular orbit, we can say that the 



number of points seen in transit is directly proportional to the 
fraction of time that the system spends in eclipse, and given a 
total number of observed data points Nor, Nr is thus 



Nr=No 



TTfl 



(14) 



where b is the impact parameter of the transit in units of the 
stellar radius, b = (a//?*) cos/, / is the orbital inclination, and 
we have assumed that R^, « a. We note that this does not 
hold for correlated noise, see I Ai grain & Pon^ (120071) . 

Taking equation[T4]and putting it back into our original for- 
mulation for x^ gives 

2 

(15) 
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R* 
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We define a new variable, xL^ as the value of equation ( fTsl l 
when ^7 = 0, corresponding to the x^ of an equatorial transit 
about a given star This gives us a new form that makes it eas- 
ier to determine when a system will have a transit geometry 
sufficient to be detected by the given observing set-up: 

x',VT^^>xL- (16) 

Rearranging, we find that we are able to detect transits up to 



„(M,.,p)=^/l-(^ 

/ieq 



(17) 



Therefore, bmax is a function of M, r, and /? as a result of 
the relations contained within the xL term. 

Finally, we can now deduce an explicit expression for Ps/n- 
Given that detectable transits can only occur between < ^ < 
bmax, we may simply integrate the normalized probability den- 
sity function of allowable impact parameters from to b„,ax to 
find a value for Ps/n- 

Given that b = (a/R^,)cosi, the maximum value of the im- 
pact parameter for a circular orbit (different from the max- 
imum at which we will detect transits) is b = a/R^,. The 
normalization for the impact parameter's probability density 
function is therefore 



1: 



Pt(b)db. 



(18) 



Since impact parameters are uniformly distributed^, this gives 
the solution Pb = Rt,/a. It should be noted that this form of 
Pb is only valid in the case of circular orbits. For the period 
ranges probed so far by transit surveys, this is not a very re- 
strictive assumption, since only ^ 10% of the known tran- 
siting planets have eccentricities inconsistent with zero. At 
longer peri ods, however, the effects o f eccentricity must be 
considered ( lBarnesll2007t lBurkell2008l) . 

If we now integrate this over the range of impact parameters 
that yield detectable transits, we find 



Ps/N(M,r,Rp,p)= 




(19) 



^ The distribution of cos i is uniform, wliicli follows from the assumption 
that planetary systems are randomly oriented. The distribution of cos i may 
be derived by considering the normalized angular momentum vector of the 
planetary orbit, which is free to point in any direction, and thus has a param- 
eter space that forms an imaginary sphere around the parent star. Particular 
values of cos i are "hoops" along the surface of this sphere, and the ratio of 
their area to that of the total sphere is the probability that cos i will take that 
value. 
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the Kepler survey, where we also examine the detectability of 
transiting planets over a range of radii and orbital periods. In 
both cases, we consider a fixed planet radius, so that 
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Fig. 1. — Comparing the analytic and exact formulations of the window 
probability. The top panel shows the window probability for detecting 2 tran- 
sits over a 30 day observing run with no nights lost and 8 hours of observing 
per night. The bottom panel uses the same situation, but is instead for a min- 
imum of three transits. 



2.1.2. />,,,„ 

The window probability (or window function) for a transit 
survey describes the probability that a transiting planet with 
a certain period will be visible to the survey, given the obser- 
vational cadence of the survey (i.e. - X hours of observing 
per night for Y nights). While many numerical routines exist 
that can exactly calculate the window probability for a given 
observing set-up, for our purposes calculating an exact nu- 
merical answer is generally too time intensive to be practical. 

We therefore developed an analytical framework that allows 
for the approximate calculation of a window probability. This 
has the virtue of requiring much less computing time. This is 
gained at the cost of a small loss of accuracy. 

We give a full description of our analytic window proba- 
bility in Appendix B. Briefly, we assume randomly sampled 
observations, which leads to simple expressions for the aver- 
age window probability. This implicitly ignores the effects 
of aliasing at integer and half-integer day periods. Thus we 
are not able to address issues related to non-uniform or biased 
a posteriori distributions o f planet parameters resulting from 
the effect of aliasing (e.g. iPont et a l. (2006)). Nevertheless, 
the analytic solution does follow the "outline" of the exact 
calculations, and closely tracks the moving average value of 
the exact calculations sufficiently accurately for our purposes, 
as is shown in Figure 1 . 



2.2. The Frequency of Planets 



df(R„^p) 
dRp dp 



df(Rp,p)/{dRp dp) denotes the probability that a given star 
has a planet with radius Rj, and orbital period p. In the exam- 
ple surveys that we consider later, we chose to focus on the de- 
tection of short period Very Hot- and Hot-Jupiters (VHJs and 
HJs) with orbital periods of 1 to 5 days, as this is the regime 
being probed by these surveys. The exception of this focus is 



df(Rp,p) 
dRp dp 



:kip)f(p)S(Rp-R'), 



(20) 



where k(p) is the correct normalization for a given period 
range, f(p) is the distribution within this range, and R' is the 
planet radius. For VHJs and HJs we use R' = lARj^p, the 

average radius of the transiting planets found to date^ . In gen- 
eral, we have left the exact distribution of /(/?) to be assigned 
as is deemed appropriate. For the VHJ and HJ period distribu- 
tions described above, we adopt a distribution that is uniform 
in period, f{p) = 1 . 

The frequency of Hot Jupiters (P < 5 days, Mp = 0.1 - 
10 Mjup) is commonly t aken to be ^ 1 % fro m the results of the 
RV planet surveys (e.g. jMarcy et al.l2005 a). However, adopt- 
ing frequencies from RV surveys can lead to biased estimates 
of the expected yields from transit surveys, because RV sur- 
veys constrain the frequency of planets as a funct ion of mass, 
rather than radius. Furthermore, as argued by iGould et alJ 
(2006b), most RV surveys are more metallicity biased than 
transit surveys, and therefore the planet frequency found by 
RV surveys will be correspondingly higher, since planet oc- 
curre n ce is an increasing function of metallicity (.Santos et alj 
I2004t iFischer & Valentil [2005). RV surveys are metallic- 
ity biased because the target samples for these surveys are 
magnitude-limited. Since metal rich stars are brighter at fixed 
color, they are overrepresented in magnitu de-limited samples 
with respect to a volume-limited sa mple (IFischer & Valentil 
120051) . Indeed, iMarcv et al] (l2005bl) note fliat because they 
use a magnitude-limited target list, their target stars have a 
definite metallicity bias. Transit surveys, on the other hand, 
are typically signal-to-noise ratio limited. At fixed color, 
metal-rich stars are both brighter and bigger These two ef- 
fects on the signal-to-noise ratio roughly cancel, and as a 
result transit surveys are very weakly biased with respect to 
metallicity (see Section 8.1 of Gould et al. 2006b for a quan- 
titative discussion.) 

Indeed, there is evidence for the effects of this metallicity 
bias in the comparison of various surveys. From the Cali- 
fornia and Carnegie Keck magnitud e-limited and met allicity- 
biased (Fischer & Valenti 2005; M arcv et al] l2005bl) survey 
sample, Cummi ng et al.l ([2008) find a frequency of '--' 1.3% 
for planets with P <5 days a nd M > O.lMjup- O n the other 
hand, f or the volume-limited ( TUdrv et aLlHoOObh CORALIE 
sample. lUdrv & SantosI (|2007|) report a frequency of ~ 0.8% 
for M > O.lMjup and f < 1 1 .5 days (a < 0. 1 AU). Assuming 
that ^ 75% of these planets have P < 5 days, this corresponds 
to a frequency of 0.6% for P < 5 days and M > 0.2My„p, a 
factor of ^ 2 lower than that found for the Keck RV survey. 
iGould et ail ( l2006bl) examined the OGLE-III transit survey, 
and found a frequency of ~ 0.5% for P <5 da ys, which was 
subse quently confirmed and strengthened by (F ressin et alj 
120071) . Thus, as would be expected based on the argu- 
ments above, the frequencies inferred from transit surveys and 
volume-limited RV surveys are similar (and entirely consis- 
tent considering Poisson uncertainties and uncertainties in the 
mass-radius relation of giant planets), whereas they are sub- 
stantially lower than the frequencies inferred from magnitude- 
limited (and metallicity-biased) RV surveys. 

' http://www.inscience.ch/transits/ 
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For the VHJ and HJ predictions that we later make us- 
ing T rES, XO, an d Kepler , we will adopt the normalizations 
from i Gould et alj (l2006bh of k(p) = 1/690 for VHJs with 
orbital periods uniformly distributed between 1-3 days, and 
k{p) = 1/310 for HJs uniformly distributed between 3-5 day 
periods. By using statistics from OGLE-III, instead of the RV 
statistics, we are avoiding the difficulties discussed above, and 
furthermore any unaccounted-for biases inherent in photomet- 
ric surveys will be preserved, and therefore any systematic ef- 
fects unknown to us will still be accurately captured in the 
final predictions. In other words, it is an "apples-to-apples" 
comparison, instead of "apples-to-oranges" using RV data. 

2.3. The Present Day Mass Function 
The next step is to integrate equation ^ over mass: 

-—J^=p4r,l,bycosb p,^„(M,r)—-dM (21) 

dr dl do Jm„„ dM 

Having previously derived an expression for Pdei, we need 
only determine the appropriate formulation for the PDMF, 
^. Using the results of iReid et all (|2002|) . who used the 
Hipparcos data sets along with their own data from the Palo- 
mar/Michigan State University survey to create a volume- 
limited survey out to 25 pc, we adopt 



dn _ f M^ 

dM \Mq 

dn _ /M, 

dM \Mq 



-1.35 



-5.2 



for 0. 1 <M,/Mq< 1, (22) 



for 1 <M,/Mq, 



where knorm is the appropriate normalization constant. A fuller 
discussion for our decision to use this particular form for the 
PDMF is in Appendix C. 

Adopting a stell ar mass den sity in the Solar Neighborhood 
of 0.032Mq pc"^ (iReid et al.ll2b02i) . we find a normalization 
of k„orm = 0.02124 pc"^. 

2.4. Galactic Structure 

While our normalization of the mass function properly de- 
scribes the density of stars in the Solar Neighborhood, be- 
yond about 100 pc we must account for variations in stellar 
density due to Galactic structure. Additionally, as mentioned 
previously, we must also take into consideration the effects of 
extinction due to interstellar dust. This is the final stage of 
integration that needs to be done on equation ([T]i: 



Nde, = 



Ji, 



I max /•*»,!, 



Pdet{r)p*{r,l^b)r cosbdrdldb. 



(23) 



2.4.1. Galactic Density Model 

We adopt the Bahcall & Soneira model for the Galactic thin 
disk ( Bahcall & Soneira .1980). The model treats stellar den- 
sity in the thin disk as a double exponential function involving 
both the distance from the Galactic center and height above 
the Galactic plane. Specifically, the mass density relative to 
the local density is given by 



p^,(r,l,b) = exp 



d — dg, 
hd,* 



K 



(24) 



where dgc = 8 kpc is the distance from the Galactic Center. 
hd^f is the scale length, and /z,^ is the scale height. 



We adopt a value of hj,, = 2.5 kpc for the scale length of 
the disk. For the vertical scale height, we use the following 
relation, wh ich is a slightly modified form of that originally 
proposed bv lBahcall & Soneiral ( Il980h : 



/i,-* = 90 



forMv <2, 



= 90-1-200(^^^^1 for2<Mv<5, 



= 290 



for 5<Mv 



We calculate the absolute V-band magnitudes of stars of a 
given mass using the blackbody assumptions, and the mass- 
radius and mass-luminosity relations described in Appendix 
A. 

We defer the complete discussion of our reasoning behind 
our choices for hd.* and /i,^ to Appendix D. The exact values 
for these parameters are not well constrained in the literature, 
and we show the effects of using different values in the Ap- 
pendix. We note that within the range of values found in the 
literature, the final prediction results may differ by as much as 
20%. 

We only treat the thin disk in our modeling because we de- 
cided that an expanded analysis that also looked at stars in the 
thick disk, halo, and bulge would not add a significant number 
of detections to our final result. Both the thick disk and the 
halo have much lower volume densities than the thin disk at 
distances typically surveyed by transit searches. In the same 
vein, the bulge is too far away. These populations are also 
typically older than the thin disk, and so have a fainter turn- 
off point, meaning that these stars are more difficult to detect. 
That being said, for deeper transit surveys, such as some of 
the all-sky surveys, these populations may play an interesting 
role. In that case, it would be fairly simple to add them into 
our model. 

We chose the Bahcall model instead of other more detailed 
models due to speed and portability considerations. While 
others, such as the Besancon models jRobin & Crezel [l986). 
take into account features such as density variations from the 
spiral structure of the Galaxy and the Galactic Bulge, they 
typically require over 30 minutes to simulate a given star field. 
Additionally, many are only accessible through webpage in- 
terfaces, and therefore difficult to include in a self-contained 
prediction algorithm. 

We do note that comparisons between our Galactic model 
and the Besancon model revealed star count differences of 
only 10% in the Kepler and TrES fields over a range of mag- 
nitude limits. 

2.4.2. Interstellar Extinction Modeling 

We also model the density of interstellar dust as a double 
exponential 



Pdust{r,l,b) = exp 



d-d. 



gc 



hd,c 



I. dust 



(25) 



Which gives the dust density at a distance r along a particular 
line of sight {l,b) relative to the density in the Solar Neigh- 
borhood (i.e. - z = Q, d = dgc). To calculate the extinction in 
a given bandpass for a star at a given distance, we then inte- 
grate equation (1251 ) over r to arrive at the effective path length 
through an equivalent Solar Neighborhood dust density. 



Qu.it= / Pdust(r,l,b)dr 
Jo 



(26) 
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The V-band extinction in magnitudes, Ay, is given by this 
path length multiplied by the appropriate normalizati on con- 
stant: Ay = kdust,\idust- We adopt kdusi,v = 1 mag kpc"' (IStruvd 
1 18471) . We determine the extincti on in other bands by using 
the extinction ratios (Ax/ Ay) from lCoxl (l2000l) . assuming that 
Rv = 3A. 

3. ASSUMPTIONS 

It is important to keep in mind that in the preceding deriva- 
tions, we made several assumptions about the nature of ex- 
trasolar planetary systems, as well as some simplifications to 
keep the math and computing time more manageable. 

1 . Detectability is Best Measured Through ~x^: While S/N 
is a better metric than photometric precision to determine if 
a transit will be detectable, and the y^ parameter is a good 
proxy for S/N, actual transit surveys use a wide variety of tests 
to detect transits. Many tests, such as the visual evaluation of 
the data by a human being, are difficult to simulate numeri- 
cally. Thus, using just a y^ test is a simplification, as human 
inspection may impose a much higher effective S/N threshold. 
As we show for typical surveys, the number of detections is 
a strong function of the S/N threshold. Thus a proper deter- 
mination of the effective S/N threshold imposed by all cuts is 
necessary before any robust conclusions can be drawn about 
the rate of detections relative to expectations. This is a funda- 
mental - and unavoidable - limitation of our study. 

2. No Stellar Limb-Darkening: We have neglected to 
include terms for stellar limb-darkening because, as stated 
above, to first-order its effects are canceled by th e the effects 
of ingress and egress on the transit lightcurve (iGould et al.l 
|2006b). 

3. Stars are Blackbodies: We assume that the radiated en- 
ergy of a star can be described as a blackbody. While in 
general this is a fair approximation (to within 0.5 mag, see 
Figures 11 and 12), it does not hold true for the bluer band- 
passes we consider (U-band in particular) and for very low- 
mass stars. 

4. Simplified Galactic Structure: While we do treat varia- 
tions in stellar density and interstellar dust densities in a gen- 
eral sense, our formalism is a simplified double exponential 
model that does not account for the spiral structure of the 
Galaxy or patchiness in the dust distribution. Therefore de- 
tailed, accurate predictions for an individual field may be bet- 
ter serv ed by incorporating m ore sophisticated Galactic mod- 
els (e.g. iFressin et all (l2007h ). 

5. No Binary Systems: While we only consider main se- 
quence stars, we do not model the statistics of binary systems; 
in our formulation all stars are treated as single. Though this 
is clearly not the case, it was outside the scope of this project 
to go deeply into binary frequencies, the ability of planets to 
form in such systems, and whether or not such planets would 
be detectable. 

However, we can make a rough estimate of the expected 
magnitude of the effect of binaries. Firs t, we assume that 1/3 
of ou r stars have no binary companion (.Duquennoy & Mayoj 
1199 Ih . We assume the remaining 2/3 of the stars are in bi- 
nary systems. A typical wide-field survey has a PSF size of 
^ 20", which corresponds to a separation of ^ 4000 AU, and 
a period of log(f /day) ^ 7.8 for an equal-mass solar-mass bi- 
nary. Adopting the p e riod d istribution for binary stars from 
iDuquennov & Mavon (Il991h . we estimate that -- 90% of the 



binary systems are unresolved. Of these unresolved binaries, 
we conservatively assume that those with separations < 5 AU, 
or log(f /day) ^ 3.5, do not host hot Jupiters because the bi- 
nary components are too close. 

Thus we estimate that ^ 62% of unresolved binaries can 
host hot Jupiters. To estimate the detectability of transiting 
planets in these, consider a binary where F\, R\, M\ and F2, 
R2, M2 are the flux, radius, and mass of the primary ("1") and 
secondary ("2"), respectively. Additionally, define i = F2/F1, 
TZ = R2/R1, and q = M2/M1 < 1. Then, for fixed planet ra- 
dius and semimajor axis, assuming photon-limited uncorre- 
cted uncertainties, the signal-to-noise ratios if the stars are re- 
solved are (S/N)i ex R^'^fI'^ and (S/N)2 ex R'^'^F^'^, while 
if the stars are unresolved the signal-to-noise ratios for transits 

across each component are (S/N)b_i cx /?[ Fj (1 +^)"'/^ 

and (S/N)b^2 oc R'^'-F^'^il +tr^l^. Assuming a constant 
volume density of stars and no interstellar ex tinction, the 
numb er of detected planets scales as (S/N)"^;^^ dGould et al.l 
l2003h . Thus the number of planet detections for a popu- 
lation of binaries in the unblended case can be written as 
Nq cx (S/N)j +/(^)(S/N)|, where f{q) is the mass ratio distri- 
bution. Similarly, for the blended case, this is A^^ oc (S/N)^ j + 
/(^)(S/N)| 2- Using the scalings for the signal-to-noise ratios 
just derived, we evaluate the ratio of the number of detections 
in the blended case, relative to the number of detections one 
would expect if they were unblended. 



Nn 



= {\+£) 



-3/2 



\+f{q)e"R-'^l^ 



(27) 



A^o ^' ' ' i+f{q)Pi^n-'^i^' 

We adopt i = (f , which is roughly appropriate for the R- 
band (see Fig. 11), and 7^ = q'"-^ (Eq. IA4I) . We consider 
two different distributions f or f(q ^. First, we adopt f{q) 
from iDuquennoy & Mayod (1199 Ih . which is weighted to- 
ward low mass ratio systems. Second, we consider a distri- 
bution with a substantial population of 'twins' with J = 1- 
Specifically, we choose f(q) cx 1 + 3.5exp[-0.5(^- 1) /cr^], 
with a^ = 0.2, w h ich ap proximates the distribution found by 
iHalbwachs et alj (l2003h . WTien we then average Equation 
dZTl l over q, weighting by f(q), we find (Nb/N)) ~ 0.89 for 
the DuQuennov & Mavo r (19911) mass ratio distr ibution, and 
(Nb/Nq) ~ 0.68 for the IHalbwachs et all (l2003h mass ratio 
distribution. Thus, we estimate that roughly 70-90% of plan- 
ets in unresolved binaries will still be detectable. The basic 
point is that, because the mass-luminosity relation is so steep, 
only those binaries that have ^ ~ 1 result in a substantial re- 
duction in the signal-to-noise ratio. Low mass ratio com- 
panions have little affect on the abilit y to detect transit sig- 
nals. F or the mass ratio distribution of IDuquennov & Mavoj 
( Il991h . these provide the bulk of the companion population, 
and so in this case the effects of binarity are suppressed even 
further. 

Finally, we can summarize these various components to 
provide a rough estimate of the fraction of planets that will 
be detected when we consider binary systems, relative to 
our fiducial assumption that all stars are single. We find 
1 /3 + (0.1 + 0.62 (A^b/M)))x (2/3) -0.68 -0.77. Thus we ex- 
pect that our estimated yields are systematically high by at 
most 30%. This is generally of order or smaller than other 
sources of uncertainty in our predicted yields. 

We note that a corollary of our assertion that binaries have a 
small effect on the yields of transit surveys is that a substantial 
fraction of the planets that have been found in wide-field tran- 
sit surveys should be in binary systems. Only two planets in 
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TABLE 1 
TrES Survey Parameters 


Band 


550-700 nm 


mn Limit 


pal3 


Observation Time 


Varies 


Exposure Time 


90 seconds 


CCD Read Time 


25 seconds 


Telescope Diameter 


0.1m 


Tlirougliput 


0.6 


PSF (FWHM) 


25" 


Opix 


10" 


CCD Full-Well Depth 


1.2xl05e- 


Field of View 


36 deg^ 


Galactic Longitude 


Varies 


Galactic Latitude 


Varies 



TABLE 2 

Properties of the Simulated TrES Survey 

Fields 



Field 


Observation Time 


Galactic Coordinates 


AndO 


765 hrs 


[126.11,-015.52] 


Cygi 


609 hrs 


[084.49, -1-010.28] 


CasO 


628 hrs 


[120.88,-013.47] 


Perl 


643 hrs 


[156.37, -014.04] 


UMaO 


497 hrs 


[168.87, -1-047.70] 


CrbO 


537 hrs 


[053.49, +048.92] 


Lyrl 


582 hrs 


[077.15, +017.86] 


Andl 


435 hrs 


[109.03, -017.62] 


And2 


332 hrs 


[115.52,-016.21] 


TauO 


406 hrs 


[169.83,-015.94] 


UMal 


500 hrs 


[156.32, +054.01] 


Herl 


592 hrs 


[063.26, +026.42] 


LacO 


686 hrs 


[097.79, -008.86] 



confirme d common proper-motion binary syst ems have been 
reported dBurke et al.ll2007t : iBakos et al]|2007 ). however sys- 
tematic surveys for stellar companions to the host stars of the 
known transiting planets have not been performed, and there- 
fore the true incidence of binary systems amongst the sample 
of transiting planets is not known. 

4. RESULTS 

Using the preceding formalism, we simulated four different 
surveys modes used by transit searches: point-and-stare, drift- 
scanning, space-based, and all-sky surveys. For the point- 
and-stare and drift-scanning cases, we use the TrES and XO 
transit surveys as our examples, respectively. Both surveys 
have been successful in their hunt for planets, and so we are 
able to directly compare the results of our calculations with 
the actual number of detections. For VHJs and HJs, we sim- 
ulate the detection of pl anets with R „ = 1 . IRjup, and with the 



period distribution from lGould et al.l (12006b.) that we describe 
in Section 2.2. 

We demonstrate our method for space-based surveys by 
simulating the upcoming Kepler mission, and comparing our 
predictions with other predictions for Kepler in the literature. 
Similarly for the all-sky mode, we examine the future planned 
observing programs of LSST, SDSS-II, and Pan-STARRS and 
make predictions for each. 

For each mode, we have tried to approximate the ef- 
fect of red noise on detectability by increasing our S/N cut 
dPont et al.li2006.) . To show the effects of a more accurate red 
noise calculation, in Section 4.5 we evaluate its effects on the 
TrES survey. 

4.1. Point-and-stare - TrES 
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Fig. 2. — The distribution of predicted TrES detections in the Lyrl field 
with and without red noise, and assuming a brighter magnitude limit than our 
fiducial value oi mr < 13. The location of the stars TrES-1 and -2 are also 
shown. TrES-2 Hes within Lyrl, while TrES-1 is in the neai'by LyrO field. 



Point-and-stare observing is the traditional observing mode 
of tracking one field of sky continuously for all, or part, of 
the night. There are several operational photometric transit 
surveys that use this particular method of obs erving, includ- 
ing H ATnet ('Bakos et al."2004). Super Wasp (iPollacco et alj 
120061) . and KELT (Pepper et al. 2007), but we chose to model 
is the Trans-Atlantic Exoplanet Survey (TrES). TrES is com- 
posed of a network of three telescopes: STARE on the Ca- 
nary Islands, PSST at Lowell Observatory in Arizona, and 
Sleuth on Mt. Palomar, California. All are small, 10 cm 
aperture, wide-field (6°), CCD cameras that operate in uni- 
son to observe the same field of sky nearly continuously 
over one to two month peri ods. The surv e y spec ifics are 
described in more detail by iDunham et al.l ( |2004|) and by 
iBrown & Charbonneaul (|1999'). The relevant quantities for 
our modeling are listed in Table 1 . 

We fit to the RMS plots posted on the Sleuth observing web- 
site** to determine the background sky brightness in the TrES 
bandpass (550-700nm), using the TrES team's stated photon 
throughput of 0.6, and we estimate that the TrES bandpass has 
a sky background of 19.6 mag arcsec"^. 

In our simulations, we treat the three TrES telescopes as 
one single instrument. This is clearly an oversimplification. 
While these telescopes are quite similar, there are differences 
between the data quality and quantity obtained from each site. 
For example. Sleuth must contend with a higher sky back- 
ground due to its proximity to urban areas, whereas the frac- 
tion of clear nights for PSST is typically lower than for the 
other two sites. However, accounting for these differences 
is difficult, and furthermore their effect on our predicted de- 
tections is small compared to the other uncertainties in our 
model. We therefore ignore these effects. In dealing with 
the window probability of the TrES network, we were greatly 
aided by the plots posted on the Sleuth website which showed 

^ http://www.astro.caltech.edu/ ftod/tres/sleuthObs.html 



Predicting the Yields of Photometric Surveys 



TABLE 3 
Predicted TrES Detections for mr < 13, S/N> 30, and Rp -- 



I.IR, 



Jup 





AndO 


cygi 


CasO 


Perl 


UMaO 


CrBO 


Lyrl 


Andl 


And2 


TauO 


UMal 


Herl 


LacO 


Total 


VHJs 
HJs 


0.38 
0.34 


0.37 
0.30 


0.34 
0.29 


0.34 
0.28 


0.18 
0.16 


0.21 
0.18 


0.32 
0.27 


0.24 
0.18 


0.19 
0.09 


0.22 
0.15 


0.17 
0.16 


0.21 
0.19 


0.41 
0.34 


3.59 
2.91 


Both 


0.72 


0.67 


0.63 


0.62 


0.34 


0.39 


0.59 


0.42 


0.28 


0.37 


0.33 


0.40 


0.75 


6.50 



TABLE 4 
Predicted TrES Detections 

FOR Mr < 13 AND Rp = l.lRjup 





VHJ 


HJ 


Both 


S/N> 20 


6.06 


5.72 


11.78 


S/N>25 


4.68 


4.05 


8.73 


S/N>30 


3.59 


2.91 


6.50 


S/N>35 


2.77 


2.13 


4.90 


S/N> 40 


2.16 


1.59 


3.75 



TABLE 6 

Predicted TrES Detections for 

S/N> 30 AND MR < 13 





VHJ 


HJ 


Both 


Rp = 0.9Rj,p 


1.66 


1.17 


2.83 


Rp = 1.0Rjup 


2.52 


1.90 


4.42 


Rp = 1-lR/ap 


3.59 


2.91 


6.50 


Rp = i.2Rjup 


4.53 


3.88 


8.41 


Rp = \..3Rjup 


5.52 


5.04 


10.56 



TABLE 5 
Predicted TrES Detections 
for S/N> 30 and Rp = l.lRjup 



TABLE? 
XO Survey Parameters 





VHJ 


HJ 


Both 


MR < 12.0 


2.13 


1.95 


4.08 


niR < 12.5 


2.95 


2.58 


5.53 


niR < 13.0 


3.59 


2.91 


6.50 


niR < 13.5 


3.86 


3.01 


6.87 


niR < 14.0 


3.94 


3.03 


6.97 


niR < 14.5 


3.96 


3.05 


7.01 



the window probabiHty for each field, as well as the total net- 
work observing time each field had received. 

We simulated thirteen different TrES fields using this 
methodology. The network observation times and the loca- 
tion of each field are shown in Table 2. The data for this table 
was collected off from the Sleuth observing website. The ex- 
pected detections for each of the fields (and the total number 
of detections over all simulated fields) are shown in Table 3. 
The effects of our Galactic structure model can clearly be seen 
in distribution of expected detections. Both the UMaO and 
CrbO fields are substantially above the plane of the Galaxy 
(b « 48°), and thus show a much lower detection rate than 
those fields directed towards denser regions of the sky. 

Figure 2 shows the distributions of detections over various 
parameters in the Lyrl field. For Figure 2, we show two distri- 
butions, one assumes white noise in the photometry, the other 
assumes 3 mmag of red noise. For the white noise estimates, 
we used a signal-to-noise limit of S/N > 30 and a magni- 
tude limit of niR < 13. Lyrl is the field in which TrES -2 was 
discovered, and is adjacent to LyrO, which contains TrES-1. 
We therefore have plotted the location of these two stars on 
the distributions shown in Figure 2. Both stars reside close 
to the peaks of our parent star mass and radius distributions. 
The third and fourth planets discovered by TrES, TrES-3 and 
TrES -4, are not plotted in Figure 2, as they reside in the Her- 
cules fields, which are at a higher Galactic latitude than Lyra. 
Also of note is that with our assumptions half of the expected 
detections the F-dwarf range. Among the false positives en- 
countered by the TrES team, the most frequent are eclipsing 
F/M-dwarf binaries. Figure 2 offers a partial explanation for 
this, showing that roughly half of TrES' planet-like signals 
come from F-dwarf s. 

Our predicted number of detections for the selected TrES 
fields of Ndei =6.50 covers only half of the fields observed by 
TrES, so we may roughly extrapolate this result to 13.00 de- 



Band 


400-700 nm 


mv Limit 


12 


Observation Time 


« 400 hrs 


Exposure Time 


54 seconds 


CCD Read Time 


n/a 


Telescope Diameter 


0.11 m (two) 


Throughput 


0.16 


PSF (FWHM) 


45" 


0p,. 


25.4" 


CCD Full-Well Depth 


1.5x10^6- 


Field Area 


7.2 X 63 deg^ 


Galactic Longitude 


Varies 


Galactic Latitude 


Varies 



tections over all of the fields. This is roughly three times the 
four planets that TrES has actually detected, though we hesi- 
tate to draw strong conclusions from this comparison. While 
we have the luxury in our predictions of imposing draconian 
signal-to-noise ratio and magnitude cuts in the data, in real 
life the TrES detections are identified in data using more flex- 
ible selection effects. A better comparison between our pre- 
dictions and the actual results of TrES (or any survey) would 
therefore require a better understanding of the selection ef- 
fects in the observed data set. 

To see how different selection criteria change our predic- 
tions for TrES, Tables 4-6 show the effects of changing the 
S/N limit, the magnitude limit, and the radius of the targeted 
planets on the detection numbers. 

A changing magnitude limit is particularly interesting, 
since all four of the TrES detections have been around stars 
with mfi < 12. This might indicate that the effective magni- 
tude limit of the TrES survey is brighter than what is described 
in the published literature. If this is the case, then using a limit 
of OT« < 12 lowers the number of expected detections in our 
simulated fields to N^et = 4.08, and total detections to 8.16. 

4.2. Drift-scanning - XO 

The XO survey is a wide-angle ground-based photometric 
survey that uses drift-scanning for its observations. The XO 
telescope is a pair of 1 1 cm cameras in place at the Haleakala 
Observatory in Hawaii. It monitors thousands of stars with 
my < 12 in six fields spaced evenly around the sky in right 
ascension^. The cameras scan within these fields, sweeping 

' 00 hr, 04 hr, 08 hr, 12 hr, 16 hr, and 20 hr. 
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TABLE 8 
Predicted XO Detections eor my < 12, S/N> 30, 

AND ftp = l.lRjup 



TABLE 9 
Predicted XO Detections 

FOR my < 12 AND Rp = l.iRj„„ 



00 hrs 


04 hrs 


08 hrs 


12 hrs 


16 hrs 


20 hrs 


Total 


0.34 


0.35 


0.34 


0.30 


0.32 


0.37 


2.02 



in declination from 0° to +63° every ten minutes. Typically, 
this provides over 3000 observations per star in each of the 
XO field s. More on the specifics of the XO system can be 
found in iMcCullough et a l. (2005). The relevant parameters 
for our predictions are listed in Table 7. 

To determine the throughput of XO and the background 
sky magnitude in the XO band (4 00-700 nm), we fi t to th e 
RMS diagram shown in Figure 8 of McCuUough et al] (l2005h . 
With the addition of their stated calibration and scintillation 
error, we were able to exactly reproduce the relations shown 
in that figure using a background sky magnitude of 19.4 mag 
arcsec"' and a system throughput of 0. 16. 

Additionally, we checked to ensure that the values we 
used for the observation time and the number of nights ob- 
served generated an analytic win dow prob ability similar to 
that shown in Figure 7 of McCullo u^h et all {2005). Our abil- 
ity to do this was greatly aided by the XO team's description 
of how many nights they lost to weather and mechanical trou- 
ble. 

We simulated each of the six XO fields for our predic- 
tions. The results for each field using the base parameters 
of S/N> 30, my < 12 and Rp = l.lRjup are shown in Ta- 
ble 8. The first planet disco vered by the XO survey, XO- 
Ib, (IMcCullough et al.ll2006l) is located in in the field at 16 
hrs, and so we have also plotted the distribution of detections 
wit hin this field over a variety of par ameters in Figure 3. XO- 
2b (lBurkeetal.ll2"007h and XO-3b (iJohns-KruU et all l2007h 
are in noticeably different fields (08 hrs and 04 hrs, respec- 
tively), so were not included on the plot. 
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Fig. 3. — The distribution of predicted XO detections in the field at 16 
hrs plotted against various parameters. Three cases are shown: our fiducial 
S/N> 30, mv < 12 case, S/N> 20, my < 12, and S/N> 20, my < 11.5. Also 
shown is the actual location of the star XO- 1 , which is located within this 
field. 





VHJ 


HJ 


Both 


S/N>20 


3.45 


2.41 


5.86 


S/N>25 


2.20 


1.27 


3.47 


S/N>30 


1.37 


0.65 


2.02 


S/N>35 


0.83 


0.33 


1.16 


S/N>40 


0.50 


0.17 


0.67 



TABLE 10 
Predicted XO Detections for 
^jiip 





VHJ 


HJ 


Both 


my < 11.5 


1.09 


0.55 


1.64 


my < 12.0 


1.37 


0.65 


2.02 


my < 12.5 


1.50 


0.69 


2.19 


my < 13.0 


1.54 


0.70 


2.24 


my < 13.5 


1.55 


0.71 


2.46 



TABLE 11 

Predicted XO Detections for 

S/N> 30 AND my < 12 





VHJ 


HJ 


Both 


R„ = Q.9Rj„p 


0.29 


0.09 


0.38 


Rp = l.ORj„p 


0.69 


0.25 


0.94 


Rp = ilRjup 


1.37 


0.65 


2.02 


^p — ^-^^jLip 


2.07 


1.18 


3.25 


R„ = l.'iRjup 


2.93 


1.92 


4.85 



Tables 9-11 show the effects of changing the S/N limit, the 
magnitude limit, and the radii of the targeted planets on XO's 
predicted detection numbers. For our base case of my < 12 
and S/N> 30, we expect that XO will detect 2.02 planets. 
If we consider the three actual detections from XO, it seems 
that, similar to TrES, the actual S/N and magnitude limit of 
the XO survey may differ from what is in the literature. For 
example, all three XO parent stars are at my < 11.5 (XO-1, 
the dimmest, is at my = 11.3). To see how this affected our 
calculated distributions, we also simulated the field at 16 hrs 
assuming S/N> 20, my < 12, and S/N> 20, my < 11.5. The 
results are included on Figure 3. In these cases, the total num- 
ber of detections increases to 5.86 and 4.21, respectively. 

4.3. Space-based - Kepler 

To make predictions for a space-based mission, we chose 
NASA's Kepler mission, which is primarily designed to look 
for Earth-sized extrasolar planets. It consists of a single 
spacecraft that is currently scheduled to be launched into 
an Earth-trailing heliocentric orbit in early 2009. Onboard, 
Kepler will have 42 CCDs that will measure the light from 
100,000 main-sequence stars within the telescope's field of 
view. Ideally, Kepler will identify planets orbiting within the 
"habitable zone" of their parent stars which would be capable 
of possessing liquid water on their surfaces. 

Beyond looking for extrasolar Earths, the immense amount 
of photometry that Kepler is expected to produce over its four 
year operational lifetime will also allow for the detection of 
larger jovian worlds transiting other stars. Previous estimates 
of the number of VHJs and HJs expected to be found by Ke- 
pler have predicted the discovery of about 180 HJs (no VHJs 
had been detected at th e time, and so they were not treated) 
d Jenkins & Dovlel[2003l) . Clearly, given that there are only 46 
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TABLE 12 

Kepler Mission Survey 

Parameters 



Band 


400-850 nm 


^Kepler 


0.86 


my Limit 


?al4 


Observation Time 


4 years 


Exposure Time 


3 seconds 


CCD Read Time 


Negligible 


Telescope Diameter 


0.95 m 


Throughput 


.3 


PSF (FWHM) 


10" 


dpix 


3.98" 


CCD Full-Well Depth 


1.2 X 10* g- 


Field of View 


106 deg- 


Galactic Longitude (new) 


76.32° 


Galactic Longitude (old) 


70° 


Galactic Latitude (new) 


13.5° 


Galactic Latitude (old) 


5° 



currently known transiting planets, this would be a significant 
increase in the number of known planets 

To simulate Kepler with our methodology, we used the 
values in Table 12 to simulate Kepler observations out to a 
distance of 8 kpc, using a signal-to-noise ratio of SN> 7. 
We gathered the information in Table 12 from both the Ke- 
pler mission website'" as well as published descriptions of 
the satellite (Borucki et al. 2004). Three items are of note. 
First, our value for the photon detection efficiency of Kepler 
is based upon information from the Kepler website that a 12th 
magnitude G2 dwarf will show 7.8 x 10^ 7 hr~' . Secondly, the 
value we use for the extinction ratio in the Kepler bandpass, 
A-Kepier/Av = 0.861, is an average that we calculated by inte- 
grating the extinction curv e over the Kepler bandpass (400- 
850nm) (.Fitzpatrickl 1 1 999ii) . Finally, the Kepler website and 
much of the literature states that Kepler will use a dim mag- 
nitude cutoff of OTv < 14. In reality, there is no "hard" magni- 
tude cutoff at which stars will be excluded from observations. 
Indeed, Kepler is expected to observe promising stars that are 
much dimmer than my = 14 (D. W. Latham, private commu- 
nication). Therefore, for our more explicit calculations of the 
expected number of HJ and VHJ detections for Kepler, we 
adjusted the magnitude limit such that the number of stars in 
the field agreed with the expected number of Kepler targets 
(100,000 in the case of the new field). For other, more gen- 
eral, simulations, we used several different magnitude limits 
to determine what effects these had on Kepler's expected de- 
tections. 

While the location of the Kepler field is currently cen- 
tered on Galactic coordinates (/ = 76.32°, b = 13.5°), up un- 
til 2004 it was instead situated closer to the Galactic plane 
at (/ = 70°, b = 5°). The switch in the position of the field 
to higher Galactic latitudes was done in order to avoid sort- 
ing through the larger number of giant stars that are found in 
fields directed along the plane of the Galaxy. Much of the 
previously published literature examining the Kepler mission 
have used the coordinates of the old field. For the purpose 
of comparison, we therefore simulated both the old and new 
Kepler fields. 

4.3.1. General Characterisation 

Looking in the new Kepler field with a magnitude limit of 
niy < 14, we simulated the number of detections as a function 
of orbital period for Jupiter-, Earth-, and Mars-sized worlds, 

'" http://kepler.nasa.gov/ 
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assuming that every star possessed a planet with the given pe- 
riod and radius (Figure 4), and that red noise was negligi- 
ble. Interestingly, the difference between the number Jupiter- 
and Earth-sized planets detected is minimal for /:> < 10 days, 
which indicates that Kepler is able to detect all the planets 
within this radius range, assuming that they transit. Therefore, 
for my < 14 and 79 < 10 days, the S/N of transiting Jupiter- and 
Earth-sized planets is large enough that Kepler's S/N limit is 
iiTelevant, and the my < 14 magnitude cutoff is the limiting 
factor in the number of detections. For smaller planets the 
size of Mars, the S/N of their transits are in the neighborhood 
of the Kepler S/N limit for the brightest stars and shorter pe- 
riods, and so Mars-sized worlds are detected at a much lower 
rate. 

To more accurately characterize the magnitude- and S/N- 
limited detection regimes, we plotted the number of detec- 
tions as a function of planetary radius for orbital periods of 
1 day and 365 days (Figure 5). In addition in a magnitude 
limit of my = 14, we also plot the same relations for a limit of 
my = 20. One can clearly see the flattening of these curves as 
the radius of the target planets increase, indicating the increas- 
ing effect of the magnitude limit on the detection numbers. 
As is expected in the S/N limited regime, detections go a s 
R^'p (iPepper et al.ll2003t iGaudJ 2005: P epper & Gau di 2005), 
though at the lower end saturation of the Kepler CCDs serves 
to cut out some of these detections. 

4.3.2. Habitable Planets 

One of the major goals of the Kepler mission is the detec- 
tion of planets in a star's habitable zone that would be capable 
of supporting life. The habitable zone is defined as the dis- 
tance from the parent star at which liquid water water could 
exist on a planet's surface. We estimate the detection rate of 
habitable zone planets by assuming that all stars have a planet 
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with a semi-major axis equal to 

AU ~ V Lbo,Q ■ 



(28) 



This assumption is simple and likely optimistic. In reality, 
the habitable zone has a finite width, and habitable planets 
will have a range of radii. Thus, a better estimate would be 
to adopt a distribution of planet masses and radii, and inte- 
grate over the range of radii and periods considered habitable. 
Furthermore, the assumption of one habitable planet per star 
is probably optimistic, and so the actual number of detected 
habitable Earths will be smaller than we estimate, depending 
on the real frequency of Earth-like planets. For the sake of 
simplicity, however, we will adopt the assumptions above. 

With this assumption, we simulated the detection of hab- 
itable Earths in both the old and new Kepler fields (Figure 
6). Surprisingly, despite there being nearly twice as many 
main sequence stars in the old field, the number of detected 
habitable Earths remained nearly the same in the two fields. 
This occurred because at the magnitude relevant for Kepler, 
the increased star count in the old field was mainly due to the 
larger number of A and F main-sequence stars found closer 
to the Galactic plane; habitable planets are doubly hard to de- 
tect around these stars due to the lessened depth of transit, 
and the longer periods required for a planet to be in the habit- 
able zone. The exact distribution of stars by spectral type for 
the old and new fields is given in Table 13, and is discussed 
further below. For the 106,000 stars in the new Kepler field, 
we find that Kepler will be able to detect 79 habitable Earth 
radius planets, assuming that every star posses one. 

In general. Figure 6 provides an excellent example of the 
benefit of repositioning the Kepler field away from the Galac- 
tic plane. Namely, the number of large dwarfs and giant 
stars (which we have not treated) found in the Galactic plane 
can be avoided, reducing the false-positive rate, while the 
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Fig. 6. — The number of habitable Earths detected by Kepler in the old 

and new fields binned by the apparent V-magnitude of the parent star, as well 

as the results assuming varying amounts of red noise on the new field. Also 

shown are the star counts for both fields binned in a similar manner 

number of habitable Earths detected remains nearly the same 
(Gould et al. 2003 ). 

It is also possible that life may arise on non-Earth-sized 
worlds. A world the size of Mars might develop quite dif- 
ferently than it has in our own Solar System if it were closer 
to its parent star, and Jovian worlds in the habitable zone 
could possess large moons capable of supporting life. Fig- 
ure 7 shows the number of detections of Jupiter-, Earth-, and 
Mars-sized habitable worlds, assuming that every star has a 
planet, against a variety of parameters in the new Kepler field. 
For completeness, the magnitude limit was set to my > 20. 
One can see in the mass and radius plots an explicit demon- 
stration of why changing the number of bright, massive main- 
sequence stars in the field does not effect the number of habit- 
able planet detections. As was mentioned earlier all detections 
drop to zero beyond 1 .4 Mq . In the plot showing the number 
of detections binned against the absolute V-magnitude of the 
parent star we have additionally shown the effect of changing 
the apparent V-magnitude limit on the distribution of Earth- 
sized planets. 

4.3.3. Hot and Very Hot Jupiters 

Using the HJ and VHJ frequencies provided bv lGould et alj 
(!2006b'), we were able to make explicit predictions about the 
number of these worlds that Kepler should find. As men- 
tioned previously, iJenkins & Doyld (l2003h (hereafter JD03) 
have made their own estimates of Kepler's HJ detection abil- 
ity (as noted previously, no VHJs were known at the time, 
and so were not considered), thus providing a useful point of 
comparison. As before, we used a limit of S/N> 7, though for 
Jupiter-sized planets Kepler is magnitude, and not signal-to- 
noise, limited. 

At the time of the JD03 estimate the Kepler field was cen- 
tered on its old coordinates. For the purposes of comparison, 
we therefore examined the old Kepler field in addition to the 
new field. 

Tables 1 3 and 14 show the distribution of our predictions by 
spectral type. To facilitate comparison with JD03, we have in- 
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TABLE 13 
Number of Main-Sequence Stars in Kepler's Field of View 



"Air 



B 



K 



M 



Jenkins & Doyle: TTOH 35708 811)85 51347 15573 831" 

Calculated (old field): 4278 31370 76630 62891 18243 1264 
Calculated (new field): 338 5066 41094 45568 13351 942 



206558 
194678 
106342 



TABLE 14 
Expected Transiting HJs and VHJs in Kepler's Field of View 



K 



M 



Jenkins & Doyle: 10 32 71 54 14 T" 

Calculated (old field): 4.2 27.3 55.5 33 8 0.4 
Calculated (new field): 0.3 4.4 28.6 23.5 5.8 0.33 



All (HJs) 



181 (181) 
127 (72.8) 
62.9 (36) 



eluded not only the final prediction numbers for both old and 
new Kepler fields, but also the distribution of the calculated 
total number of stars in each field. As stated above, in both 
fields we set the magnitude limit of the simulation such that 
the total number of stars in the fields agreed with the number 
of expected Kepler targets. In the case of the new field, this is 
100,000 stars, while for the old field we matched the number 
of stars predicted by JD03. This corresponded to a magnitude 
limit of mv = 15.7 in the new field and niv = 15.9 in the old 
field. 

Looking at the detections, one can see in Table 14 that our 
final prediction of Ndei = 127 for the old Kepler field is lower 
than the JD03 prediction. This is primarily a result of the 
fact that the Gould et al. (2006b) results that we use for k{p) 
(see Section 2.2) posit a HJ and VHJ frequency approximately 
1/3 of the frequency used in JD03. An important difference, 
however, between our work and JD03 is our consideration of 
VHJs, which were not treated by JD03. This has the effect of 
"doubling" the number of planets we expect Kepler to detect 
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(as demonstrated by the final values in parentheses in Table 
14). We therefore arrive at detection rates for VHJs and HJs 
that are 2/3 that of JD03 for the old Kepler field, and 1/3 that 
of JD03 in the new Kepler field. For the old field, we predict 
54.2 VHJs and 72.8 HJs. For the new field, we predict 26.9 
VHJs, and 36.0 HJs. 

In the new Kepler field, we see that looking away from the 
plane of the Galaxy reduces the number of HJs and VHJs de- 
tected by half to Ndet = 62.9, which matches the overall drop 
in the star count in the new field. The distribution of detec- 
tions changes substantially, shifting more towards F-, G-, and 
K-dwarfs. This is the caused by our model of the Galaxy's 
vertical structure, which uses reduced scale heights for earlier 
stars. Thus, most of the difference in star counts between the 
old and new fields is caused by fewer B-, A-, and F-dwarfs, 
leading to proportionally smaller detection numbers around 
these types of stars. 

While the prediction number of detections in the new Ke- 
pler field is lower than the corresponding JD03 prediction 
would be for the same star count, this is accounted for in the 
same way as in the old field. Indeed, we find that our predic- 
tions are consistent with those of JD03, except for the differ- 
ence in our assumption about the frequency of HJs and VHJs. 
The eventual results from Kepler therefore offer an excellent 
opportunity to more robustly determine the actual prevalence 
of short period Jovian worlds. 

4.4. All-sky Surveys 

In addition to wide-field photometric surveys, several cur- 
rent and planned projects aim to repeatedly gather photometry 
over gigantic swathes of the sky. While not all of them have 
extrasolar planet detection as one of their goals, the large ar- 
eas covered by these surveys, together with their long dura- 
tions (from 3 to 10 years), offer the possibility of detecting 
hundreds of transiting VHJs and HJs. 

We examined three all-sky surveys, LSST, SDSS-II, and 
Pan-STARRS. The Lai'ge Synoptic Survey Telescope (LSST) 
will be a large 8.4 m telescope, with an effective aperture of 
6.4 m, situated in northern Chile. It is scheduled to finish con- 
struction sometime in 2014. Using a 9.6 deg^ field of view, the 
LSST team intends to image 20,000 deg^ of the sky repeatedly 
over the course of 10 years." 

SDSS-II is the second phase of the Sloan Digital Sky Sur- 
vey (SDSS) which began in June, 2005 and is scheduled to 
last for three years until June of 2008. It is composed of three 
distinct surveys, the Sloa n Legacy Survey, SEGUE, and the 
Sloan Supernova Survey dSako et alj|2005h . The Supernova 

' ' http://www.lsst.org 
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TABLE 15 

All-Sky Survey Parameters 



Survey 


LSST 


SDSS-II 


Pan-STARRS 


Pan-STARRS Wide 


my Limit (Sun-like) 

my Limit (M-dwarfs) 

Observation Time 

Telescope Diameter 

Throughput 


18.51 

23.13 

10 years 

6.5 m 

0.5 


15.56 

20.18 

37.5 days 

2.5 m 

0.5 


14.99 

19.61 

5 months 

1.8 m 

0.5 


12.49 

17.11 

5 months 

1.8 m 

0.5 


Field of View 
Area Surveyed 


9.6 deg^ 
20,000 deg2 


6.25 deg^ 
300 deg2 


7deg^ 
1200 deg- 


7deg^ 
12,000 deg' 



TABLE 16 
Transiting VHJ and HJ Densities Around Sun-like (G2) Stars 

(DEG^^) 



Mag. Limit 


Gal. Plane 


Gal. Poles 


All-Sky Average 


mv < 20 


4.052 


0.027 


0.800 


LSST (18.51) 


1.609 


0.027 


0.387 


mv < 18 


1.125 


0.026 


0.293 


mv < 16 


0.219 


0.025 


0.087 


SDSS-II (15.56) 


0.146 


0.020 


0.063 


Pan-STARRS (14.99) 


0.083 


0.016 


0.041 


mv < 14 


0.029 


0.009 


0.017 


Pan-STARRS Wide (12.49) 


0.005 


0.003 


0.004 


mv < 12 


0.003 


0.002 


0.002 



TABLE 17 
Transiting VHJ and HJ Densities Around M-dwarfs (M5) 

(DEG^^) 



Mag. Limit 


Gal. Plane 


Gal. Poles 


All-Sky Average 


LSST (23.13) 


2.7558 


0.0888 


0.7764 


SDSS-II (20.18) 


0.2473 


0.0398 


0.1136 


mv < 20 


0.2081 


0.0368 


0.0989 


Pan-STARRS (19.61) 


0.1422 


0.0304 


0.0725 


mv < 18 


0.0257 


0.0105 


0.0169 


Pan-STARRS Wide (17.11) 


0.0092 


0.0048 


0.0068 


mv < 16 


0.0047 


0.0015 


0.0028 


mv < 14 


0.00017 


0.00015 


0.00016 


mv < 12 


0.00001 


0.00001 


0.00001 



transiting VHJs and HJs they should discover, we first cal- 
culated the V-magnitude limits at which each of the surveys 
will be able to detect transiting Jovian worlds. The complete 
derivation of this limit is given in Appendix E. Table 15 shows 
the results of these calculations, as well as the other relevant 
parameters for each of the surveys. We calculated magnitude 
limits for detecting giant planets around both Sun-like (0.8 < 
M^/Mq < 1.2) stars and M-dwarfs (0.1 <M* /Me < 0.5) for 
S/N> 20. 

To translate these magnitude limits into the number of tran- 
siting planets each of the surveys should detect, we used our 
methodology to calculate the density of VHJs and HJs for var- 
ious limiting magnitudes, in numbers per square degree, along 
the Galactic plane, at the Galactic poles, and averaged around 
the sky. Table 16 shows the resultant values for a series of 
fiducial magnitude limits, as well as the planet densities at the 
magnitude limits of our modeled surveys for Sun-like stars; 
Table 17 shows transiting planet densities around M-dwarfs. 
Additionally, in Figure 8 we show the transiting planet density 
around Sun-like stars and M-dwarfs as a function of Galac- 
tic latitude. It is interesting to note the dip near the Galactic 
plane, which occurs as a result of extinction from interstellar 
dust. 

Knowing the angular size of each survey, we use these den- 
sities to predict the number of planets that will be seen. For 
SDSS-11 and Pan-STARRS, which are both pointed towards 



Survey is the only one of these three that will collect photom- 
etry capable of finding transiting planets; the other two survey 
modes are focused on non-variable objects within and outside 
the Galaxy. To execute the Supernova Survey, SDSS uses a 
2.5 m telescope on Apache Point in New Mexico with a 6.25 
deg^ field of view. Over the three year run of SDSS-11, the 
team expects the Supernova Survey to image about 300 deg^ 
of the sky near the Galactic poles every two nights for the 
three months a year the poles are visible. Given the observing 
cadence the survey expects to use, this translates into a total 
observing time of 37.5 days. 

Finally, Pan-STARRS - the Panoramic Survey Telescope 
and Rapid Response System - will be a group of f our 1.8 
m telescopes atop the summit of Haleakala, Hawaii ( Kaiseij 
12004). The first telescope, PS-1, is nearing completion. When 
fully operational, Pan-STARRS will embark on a 3 year-long 
"Medium Deep" survey that will observe 1200 square degrees 
of the sky near the Galactic poles for a total of five months.'^ 
We also looked at a notional "wide" Pan-STARRS survey that 
looked at lOx the amount of sky as the Medium Deep survey, 
but used the same amount of observing time. 

To characterise each of these surveys and find how many 

'^ The remaining balance of the 3 years will be taken up by the other Pan- 
STARRS science modes; see 

http://pan-starrs.ifa.hawaii.edu/project/reviews/PreCoDR/|.. 
documents/PSCoDD_l_4_design_reference_mission.pdf^ 
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Fig. 8. — The number of transiting planets per square degree around Sun- 
like stars and M-dwarfs for several V-magnitude limits. The left plots are 
differential, showing the value at a particular Galactic latitude. The right plots 
are cumulative, showing the total planets from zero degrees up to a particular 
galactic latitude. 
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TABLE 18 
All-Sky Number of VHJs and HJs Detected 



a term to include a covariance coefficient Cij between the /' 
and /'' measurements, 



:lli 





Sun-like (G2) 


M-dwarfs (M5) 


LSST 


7740 


15530 


SDSS-II 


6.0 


11.9 


Pan-STARRS 


19.2 


36.5 


Pan-STARRS Wide 


48.0 


81.6 



the Galactic poles, we used the polar densities, while for 
LSST and the wide Pan-STARRS, with their much larger and 
wide-ranging survey area, we used the all-sky average density. 
The final predictions for both Sun-like stars and M-dwarfs are 
shown in Table 18. 

We note that these predictions are based on a simpler model 
than what we have used in other portions of this paper. For in- 
stance, we do not include the effects of CCD saturation or sky- 
background. Nevertheless, our calculations should provide a 
good expectation for the range of the detection numbers that 
these all-sky surveys can expect. 

The large difference between the number of planets ex- 
pected to be detected by LSST and the other surveys is caused 
by two effects. First, SDSS-II and Pan-STARRS cover much 
smaller areas of the sky to shallower depths than LSST, and 
second, LSST is not primarily directed towards the Galactic 
poles, greatly increasing the density of stars for which planet 
detection is possible. Pan-STARRS Wide, meanwhile, has a 
substantially brighter magnitude limit than any of the other 
surveys, which prevents it from detecting as many planets as 
LSST, despite the similarities in the size of their survey areas. 
Pan-STARRS Wide does detect more planets than the planned 
Narrow survey, primarily because the latter "runs up" against 
the scale height of the Galaxy. 

The real difficulty, however, in managing the 7,740 or 
15,530 detections that we expect LSST to make will be prop- 
erly identifying and separating the false-positives from the ac- 
tual planetary transit events. Indeed, to date the wide-angle 
surveys have generated approximately 15 false positive can- 
didates for every actual transiting planet they have seen (D. 
W. Latham, private communication). If we extend this ratio 
to LSST, SDSS-II, Pan-STARRS and Pan-STARRS Wide, we 
would expect them to create 116,000, 90, 290, and 720 false 
positives, respectively around Sun-like stars, and 233,000, 
180, 550, and 1,200 around M-dwarfs, respectively. Working 
through this many transit candidates (especially in the case 
of LSST) is an intensive, time-consuming task. We would 
therefore expect the all-sky surveys, despite the high absolute 
number of predicted detections, to discover transiting planets 
at about the same rate as the operational wide-angle surveys, 
as the identification rate is predominantly limited by the abil- 
ity to make follow-up observations of candidates to determine 
whether or not it is a false positive. 

4.5. Effects of Red Noise 

As pointed out bv lPont et alj (l2006h . the general mathemati- 
cal description we have used for the S/N of a transit, equation 
(in, assumes uncorrected errors on the collected photomet- 
ric data points (white noise). In practice this may not be the 
case, as phenomenon such as changes in the seeing, airmass, 
telescope tracking, etc. introduce systematic trends into the 
observed data that are correlated on the timescales of transits 
(red noise). This has the effect of substantially decreasing the 
effective S/N displayed by a transit, t hus making detections 
more difficult. Specificallv. lPont et alJ (12006.) show that to ac- 
count for red noise, equation (|4|l must be modified by adding 



X 




(29) 



This is equation (4) in iPont et al] (l2006 f). using the variables 
we previously defined for our equation (|4|. The authors give 
a more manageable version on page 8 of their work, which 
relates the S/N with red noise to the S/N without red noise: 



x^ 



Xw 



(30) 



where Xw is the white noise x^ statistic we calculate in equa- 
tion ^,Nkis the number of data points in the single k''' transit 
(which on average is just A^,r/number of transits), and a,- is the 
red component of the noise in the data. 

To properly use this formulation , we must therefore h ave 
some idea of what value a,- takes. P ont et al.l (120061 1200 7^ in- 
fer values of 3 mmag of red noise for typical ground-based 
surveys, which can be reduced to approximately 2 mmag af- 
ter applying de-trending algorithms to the data. We there- 
fore took this tentative value of a,- = 0.002 and used it to re- 
simulate the TrES fields to see what effe ct red noise woul d 
have on detection numbers. As noted in iPont et all ( |2007h , 
we set a signal-to-noise ratio limit of S/Nr > 12, which is 
the S/Nr value of the detection of TrES- lb; TrES-2b has a 
higher value of S/Nr = 14. The results are shown in Table 19. 
Figure 2 shows (in addition to TrES with only white noise) 
the distribution of detections using 2 mmag of red noise and 
S/Nr > 12. 

Compared to a white S/N limit of S/N^ > 30, using cr,- = 2 
mmag and S/Nr > 12 drops the total number of expected de- 
tections by about 20%. Looking at the split between HJs and 
VHJs, the HJ population is disproportionately effected, which 
is a result of the relatively fewer transits HJs show, lower- 
ing the overall S/N of their transit signal. Similarly, Figure 2 
shows that red noise also cuts out the detections around higher 
mass stars, as these have a lower S/N overall. 

Without de-trending the data, a red noise of 3 mmag would 
cause an even further reduction in detections, down 60% as 
compared to white noise. 

We also examined the possible effect that red noise would 
have on the yield of habitable Earth-radius planets from Ke- 
pler mission. Given that Kepler is a space-based mission de- 
signed to detect transit signals with fractional depths at the 
level of ^ lO"'*, we would expect that the amount of red 
noise due to instrumental and environmental sources would 
be small, if it exists at all. However, this is generally diffi- 
cult to predict before launch. Furthermore, there may exist 
significant astrophysical red noise due to intrinsic variability. 

We therefore estimated the expected Kepler yield for a set 
of three values for 0-^, assuming that every star has an Earth- 
radius planet with a = A\](L/Lq)^/^ (Table 20). For a rel- 
atively large amount of red noise at fXr = 10"^, we find the 
number of habitable Earths-radius is drastically reduced from 
our baseline (uncorrelated noise) prediction, from 78.89 down 
to 14.37. This is not surprising, since Earth passing in front of 
a Sun-like star has a transit depth of only 8.5 • 10"^. For this 
amount of red noise, Kepler is restricted to detecting Earth- 
like planets only around K- and M-dwarfs. This shifts the ex- 
pected distribution of host stars towards dimmer magnitudes 
(Figure 6). For smaller values of cr,. = 5 • 10"^ and cTr = 10"^, 
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TABLE 19 
Predicted TrES Detections With 2 mmag Red Noise, S/N> 12 





AndO 


Cygi 


CasO 


Perl 


UMaO 


CrBO 


Lyrl 


Andl 


And2 


TauO 


UMal 


Herl 


LacO 


Total 


VHJs 
HJs 


0.34 
0.25 


0.33 
0.22 


0.31 
0.21 


0.30 
0.21 


0.17 
0.11 


0.20 
0.13 


0.30 
0.19 


0.22 
0.12 


0.17 
0.07 


0.20 
0.11 


0.16 
0.11 


0.20 

0.14 


0.36 
0.25 


3.26 
2.15 


Both 


0.59 


0.55 


0.52 


0.51 


0.28 


0.33 


0.49 


0.34 


0.24 


0.31 


0.27 


0.34 


0.61 


5.41 



TABLE 20 

Kepler and Red 

Noise 



(Jr 


Hab. Earths 


KH 


14.37 


5 • 10-5 


32.53 


10-5 


72.88 





78.89 



the number of habitable Earths detected is reduced to 32.53 
and 72.88, respectively. Both of these estimates also assume 
that every target star has an Earth-sized planet in the habitable 
zone, and therefore must be viewed as upper limits. 

Thus we conclude that Kepler's yield of habitable zones of 
habitable planets will be not be strongly affected by red noise, 
provided that the red noise is < 10"^ on the timescales of the 
transits, which is ^ 13 hrs(M/M0)''^^ for solar-type stars. If 
the red noise is found to be larger than this, it may compro- 
mise the primary Kepler science by significantly reducing the 
number of expected habitable planet detections. For example, 
if the red noise is (Jr = 10"^, and the frequency of habitable 
planets with r = T?® is ^ 20%, there is a ~ 5% chance that 
Kepler will not detect any. The effects of red noise on the 
habitable planet yield can be mitigated by including a larger 
number of K and M dwarfs in the target sample. Of course, it 
will also be possible to extrapolate the frequency of shorter- 
period or more massive planets, where the effects of red noise 
are mitigated, into the habitable zone. 

5. CONCLUSION 

In this paper we have developed a method to more accu- 
rately predict the number of transiting planets detectable by 
wide-field ground- and space-based photometric surveys. We 
have taken into account such factors as the frequency of gas 
giants around main sequence stars, the probability of transits, 
stellar density changes from Galactic structure, and the effects 
of interstellar extinction. We adopt a S/N detection threshold 
criterion to approximate the detectability of planetary transits. 

To test our model, we compared our predictions to two op- 
erating photometric surveys, TrES and XO. For both, we cal- 
culated the number of statistically expected detections, as well 
as the population characteristics of the host stars. For TrES, 
we overestimate the number of detections (thirteen versus the 
four actual detections), while for XO our simulations yield 
approximately the actual number of detections (two predicted 
versus three actual). The exact reason for the discrepancy be- 
tween our predictions and the actual TrES results could be 
one of many. As we discuss in the TrES section, from the 



properties of the known TrES planets, it may be that the sur- 
vey has a brighter magnitude cut off than the one we use in 
our fiducial case. If this were the case, and we changed the 
magnitude limit from m« < 13 to m« < 12, the number of ex- 
pected detections would drop to 8.16 from 13. Red noise also 
has the potential to substantially affect TrES; a specific esti- 
mate of the magnitude of red-noise in the TrES survey would 
be worthwhile. 

Indeed, it is generally hard to draw specific statistical con- 
clusions from the comparison of our results and those of the 
transit surveys because they typically do not adopt the strict 
detection criteria we have used. Promising planet candidates 
are often followed up even if they are beyond the stated S/N or 
magnitude limits of the survey. Understanding and quantify- 
ing how the survey teams select candidates is vital to appro- 
priately deriving the statistical properties of extrasolar plan- 
ets. Indeed, Tables 4-6 and 9-11 demonstrate that the actual 
predictions depend crucially on the specific magnitude and 
S/N threshold used. 

As more transiting planets are discovered, these statistical 
properties are increasingly becoming the frontier of research, 
shifting the focus of the field away from individual detections. 
Given that transiting planets are so far one of the only ways 
that we may infer the radius of extrasolar planets and their 
exact orbital properties, the statistical characteristics of this 
group of objects is one of the only foreseeable ways that the 
areas of planetary interiors, system dynamics, migration, and 
formation will acquire more data. 

In the future, besides using the preceding formalism to de- 
sign more efficient transit surveys, we would hope that a more 
systematic approach towards transit surveys will allow this 
model to be used to make more specific statistical compar- 
isons. As the number of known transiting planets grows, we 
would also expect that this formalism will be used to test dif- 
ferent distributions of planet frequencies, periods, and radii 
against those observational results. This will allow us to better 
understand the statistics of extrasolar planetary systems, im- 
prove our ability to find new planets, and help to understand 
the implications of the ones that have already been detected. 
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APPENDIX 
PHOTON LUMINOSITY OF A SOURCE STAR 



The photon luminosity in a bandpass of a given star of a given mass is calculated exactly by dividing the Planck radiation 
law for wavelength by the energy of a photon at that wavelength, multiplying the result by the surface area of the star, and then 
integrating over the wavelength range of the relevant bandpass: 



$a(M) : 



Ex 



4Tr^RldX, 



(Al) 
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Where B\(T^,) is the Planck Law, and E\= ^ is the energy of a photon. Note that this assumes that stellar radiation can be 

approximated as a blackbody, a crude assumption that still suffices for our purposes, as we show below. 

Because the widths of the bandpasses in the visible and near-infrared that we are considering are generally narrow with respect 

to ^j^, we can safely approximate this integral by the product of the central wavelength of each bandpass Ac and a bandpass 

width AA: 

_ B^(7VOAA 2 2 _ 87T^cRl\-'A\ 

5'a(m) = ^ 47r K^ = — -, (A2) 

£a, exp[/ic/ Ac^T; (M)] - 1 

which, once we select a bandpass, only leaves the temperature of the source star T^,{M) to be determined. 
Using the Stefan-Boltzmann Law, we may relate the temperature of the source star to its bolometric luminosity and radius: 

where /?» and L/,o/,* are the radius and the bolometric luminosity of the source star, respectively. We chose to treat both as 
functions of the mass of the source star, as this is an independent variable that allows the use of the present day mass function to 

describe the general stellar population.^ 

Using stellar data from Harmaneq (|1 988l) and iPoppeH (|1980|), we fit broken power laws for both the R^,{M) and Lhoi.*(M) 
relations. For the mass-radius relation, we fit for two breaks in the power law, while the fitting for the mass-luminosity relation 
was accomplished using only one break. The relations that we thus found are: 

R^{M) = Rq{^^^] a = 0.8 for M* < 1 , (A4) 



M, 







M 



a=1.5 fori <M, < 1.5, 
a=.54 for 1.5 <Mh<, 



/5 



L,„,,,(M) = L,„,,o( — 1 /3 = 2.5 forM, <0.4, (A5) 

/3 = 3.8 for 0.4 <M^. 

Equations ( IA4l l and ( IA5I ) are plotted against the 'Harmaned (119881) and'Poppeij (Il980h data in Figure 9. 

The real test, however, is how well these relations allow for the calculation of the effective temperature of a star at a given 
mass using ( IA3b . as this feeds directly into the determination of t he photon luminos ity i n a given band. Figure 10 shows a plot 
of our mass vs. temperature relation alo ng with da t a poin ts from iHarmaned d 19881) andlPopper (1980). Additionally, we also 
show the analytical relation provided by IHarmaned (Il988l) . One can see tha t our formulation of (IA4b and (IA5b affords a close 
approximation to both th e da ta and the relation provided by IHarmaned (Il988h . 

Substituting back into ( IA2b allows us to calculate the photon luminosity of the source star, as well as the absolute magnitude of 
the source star in various bandpasses. If we compare these calculated absolute magnitudes in different bands for various stellar 
masses against empirical magnitudes for those same masses (Figures 1 1 and 12) we can see that the two relations are similar, to 
within half a magnitude, especially in longer wavelengths. Shorter wavelength bands such as U, B, and to a certain extent V, have 
calculated magnitudes that diverge from the empirical relation below about a solar mass, which is an unfortunate effect of using 
blackbody relations to determine the photon luminosity, as in these higher energy bands stellar spectra increasingly do not look 
like black bodies. Nevertheless, the blackbody assumption is sufficiently accurate for our purposes. 

ANALYTIC WINDOW PROBABILITY 

Roughly, if we ignore the effects of aliasing, we may say that the window probability of seeing a minimum of N,r,mm transits of 
a planet with orbital period p over the course of A^„ nights may be broken up into four different regimes. 

First, if the average number of transits during the total time during which the target star is observed exceeds Nrr.min, then it is 
statistically "certain" that the transit will occur during the times of observation . 

Second, it may be that the entire length of the survey is long enough that there are at least A',,.,,,,,, transits, but that during the 
total time spent observing the average number of transits is less than that required for detection. Then one is not certain to detect 
Nir.min transits, and there is a decreasing window probability (for increasing planet period). 

Third, one may have A',,-,,,,,,, transits during an observing run, but only if the first transit occurs sufficiently early during that 
particular run. The probability of detecting the system then decreases at an even greater rate as the orbital period increases. 

Finally, the orbital period of system may be so great that there is no possible way in which one could see the minimum number 
of transits required for a detection, and Pmndow = 0. Mathematically then, 

n«Wov,..= l for 1 < ^J^"' ^ , (Bl) 

for < 1 < 



^tr.min J ^^tr.min ^^tr 
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< JV,, > y / (JV„-l)A+7;,-(JV„„„„-l)P 

^^tr.min J \ ^ 



for < 1 < 



= 



N. ■ P 

^ {N„-\)K + T„ 

for < 1 , 

{N,r.,nin-\)P 



V^tr.min ^ )^ 



where A = 1 day. < Ntr > is the average number of transits that will be seen during the time spent observing, and is given simply 
by: 

<N,r>=^. (B2) 

Ng is defined as the total number of "good" observing nights in a particular run, and may mathematically be described as Ng = 
N„-Niost, where Niost is the number of nights lost due to weather, technical problems, or other issues. T„ is the length of each 
night's observing time. Alternatively, one may simply view the quantity A^^T;, as the total amount of good observing time in a 
given run. 

To test the accuracy of our formulation, we integrated and found the average window probability for both the analytic and 
exact formulations linearly averaged over several period ranges (1-3, 3-5, and 5-10 days). These are shown as a function of total 
observing nights under different viewing scenarios in Figure 13. 

One can see that as a result of the analytic solution not taking into account aliasing, there is divergence between the two 
solutions as the analytic form approaches unity, as well as in the cases requiring a minimum of three transits. 



THE PRESENT DAY MASS FUNCTION 

A review of the literature shows that while there is agreement that the PDMF beyond IMq can be treated as a power law with 
a slope of a Ki 5.3, below a solar mass this agreement breaks down. This is because the low-mass end of the mass-function is 
unfortunately not as well constrained as it is for higher mass stars, reflecting the difficulty of observing dim M-dwarfs, as well as 
uncertainties in the mass-lumino sity function (Henrv et al. 1999). 

For instance, Reid et al.l (|2002|) . who used the Hipparcos data sets along with their own data from the Palomar/Michigan State 
University (PMSU) survey to create a volume-limted sample out to 25 pc, propose a traditional power law formulation of the 
PDMF: 

nM)=^ = kpDMFM-", (CI) 

dM 

with a = 5.2 ±0.4 for M >Mq, and a = 1 . 35 zt .2 forM <Mq, and kp^MF as the correct normahzation. 
In a recent review, meanwhile. IChabrien (l2003h argues for a log-normal function below a solar mass: 



*(M) = 



dn 



1 



dM M(lnlO) 



A exp 



-{\ogM-\ogM,Y 

2^2 



(C2) 



3+0.021 



with A = 0.158!^:;;^^, Mc = o.o79:,,oi6: 



andcr = 0.69!*]:jj^. 



: 0.019Mq1pc- 



For comparison, if we normalize both to the lScalol(ll986h normalization for 5 Gyr of {dn/dM)\MQ = kpoMF 
and plot the two functions against each other (Figure 14), we can clearly see the divergence both numerically and functionally 
of these two forms below a solar mass. Interestingly however, for most surveys this divergence makes little difference in the 
final number of predicted detections: typically Reid's relation gives only 3% fewer detections than Chabrier's does, if both are 
normalized at a solar mass. That this difference is so small comes from the fact that most wide-angle surveys are not sensitive 
to detections around stars with log(M) < -0.15. Hence, the region of greatest uncertainty in the PDMF is left unseen by these 
surveys. 

Since the diffe rence between th e two relations is small, and becaus e the data set used bv lReid et al.l (|2002|) is much more recent 
than that used by lChabrieJ (l2003h . we therefore chose to use the Rei d et all (l2002h form of flie PDMF. 

For the purpose of comparing the Reid and Chabrier PDMFs in Figure 14, we normalized both to the value given by'Scalo 
( 1986). Going forward, however, we change the normalization to the integrated mass density found observationally bv Reid et al. 
|2002Dof0.032Mopc-^ 



GALACTIC STRUCTURE 



The values for the scale length and scale height of the Galactic disk that we use in our modeling are /j^ , = 2.5 kpc and h~„ as 
a function of absolute magnitude. Specifically, we calculate the vertical scale height of the disk as: 



/!2,* = 90 



= 90 + 200 
= 290 



m 



for Ml/ <2, 

for 2 <My <5, 
for 5 <Mv, 



(Dl) 



which is a modified form of the relation found in iBahcall & Soneird (Il980h . who instead used a value of 320 pc for the scale 
height of stars with 5 < My ■ 
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The large uncertainties in the value of the scale length of the disk reflects the uncertainty in the relevant Uterature. In their 
initial paper. | Bahcall & Soneiral(ll980h use hri «, =2.5-3.0kpc. Since then various sources have propo sed an even hd.* =3.0 kp c 
jKent et al.lll991l) . have covered all the bases by using sever al values between h^i t, = 2. 0-3.0 kpc dOlling & Merrifieldll2001h . 
and have given the very exact result of /ij „ = 2.264 ± .09 kpc dPrimmel & Spergelll2001h . 

To see what effect different scale lengths have on our final detection numbers, we simulated observations of a 6° x 6° field 
along the Galactic plane at intervals of 10° in Galactic longitude, as well as fixing the Galactic longitude at certain points (0°, 
90°, and 180°) and varying the latitude of the fields. Figure 15 shows these results for scale lengths of 2, 2.5, and 3 kpc. At the 
extremes (directly towards or directly away from the Galactic core) the 2 and 3 kpc scale lengths differ by only 0.15 detections 
per field. We chose 2.5 kpc not only to split this difference , but a lso because more recent results place the scale length much 
lower than earlier 3 kpc estimates [e.g. - iDrimmel & Sp ergel (2001)1. 

In the case of the scale height of the Galactic disk, the literature divides between those who derive the scale height from optical 
star counts and those who rely on near-infrared observations. In the optical case, the scale height is presented as a function of 
absolute visual magnitude that varies from 90 pc to values between 250-320 pc. The near-infrared observations, on the other 
hand, lead to one single value of the scale height, generally around 290 pc. This happens because measurements in the infrared 
are dominated by older M-dwarfs that have dispersed in an even fashion away from the Galactic plane. 

To compare these different estimates, we have plotted the number of detections in a 6° x 6° field along I = 90° as a function of 
Galactic latitude in Figure 16. We used two of the optically-derived scale height relations, as well as a constant scale height of 
290 pc. One of the optical relations, the Bahcall model, is given by 

/!,,* = 90 forMv<2, (D2) 

= 90 + 230 (^^^Y^j for2<Mv<5, 

= 320 for 5 < My. 

The other optical model is taken from Allen's Astrophysical Quantities (AAQ) (ICoxll2000l) and is: 

/!;,* = 90 forMv,<2 (D3) 

= 90+ 160 (^^^^j for2<Mv<5, 

= 250 for 5 < My. 

Also plotted in Figure 16 are the number of detections above and below My = 5 for the three different scale height relations, as 
well as the residuals of the constant 290 pc and the AAQ relation against the Bahcall formulation. The exact choice of how the 
Galactic scale height is treated leads to rather different final predictions. Specifically, while the optical relations are similar, the 
constant 290 pc scale height from near-infrared observati ons leads to substan tially more detections around stars with My < 5. 

For our predictions, we have adopted the form of the iBahcall & Soneiral (|l980) and Allen's Astrophysical Quantities scale 
height relations. We find the gradated scale height in these models appealing, because it accounts for the fact that more luminous 
stars with shorter lifetimes should not be dispersed far away from star-forming regions in the plane of the Galaxy. The constant 
290 pc relation, because it is based on near-infrared data, reflects the sensitivity of those observations to less massive stars that 
are brighter in longer wavelengths. Nevertheless, we have used the 290 pc result as our value for the scale height of late-type 
stars with 5 < My, in contrast to the values used in the AAQ and Bahcall models. 

ALL-SKY SURVEYS' LIMITING MAGNITUDES 

As per equation dU, we have 

r 

^=S/N = (N,r)'^^-. (El) 

a 

For the all-sky surveys that we are considering, the number of in-transit data points observed is simply the number of data 
points acquired in each field multiplied by the fraction of time a planet will spend in transit during its orbit: 

S/N=(l:!!^^Y'i, (E2) 

\texp^ ^survey "^O. J (7 

where tots is the total observation time of the survey, Qfov is solid angle of the survey telescope's field of view, and Vtsun-ey is the 
total solid angle of the sky covered by the survey. 
Rearranging: 

, texp ^survey T^a J S/N 

For Poisson, source-noise dominated errors, we have 

'^ ~ '^source (F/I-) 

where Nsomce is the total number of photons from the source collected during t^xp- We can then write, 

\ "'/^ / n \ "' 

lexp ^A \ ( LJ scope \ . ^0.2(Hiv-mv,o) 



texp.O / \ Lf scope fi j 
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where ctq is the fiducial fractional photometric uncertainty achieved for a star with apparent magnitude myfl using a telescope 
with diameter Dscope.o and an exposure tim e of f„p,o- We have defined e\ as the efficiency of the observing setup relative to the 
fiducial case. If we rearrange equation (IE5b to 



niv = 5 log 



^exp ^X 



1/2 



D 



scope 



+ mvfi 



texp.O ) \DscopefiJ <JQ 

we may now substitute in equation ( IE4b for a. Doing so, and simplifying, we find that 



my = 5 log 



tabs e\ ©FOV R* 



1/2 



^ scope ^ 



texp.O ^survey Tl'fl / D scope fi <J0 S/N 

Which is the limiting magnitude around which an all-sky survey will detect planets. 



+ '«v,o 



(E6) 



(E7) 
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Fig . 9 . — The mass-radius and mass-luminosity relations. The data points from lHarmaned <1988l) and lPoppe3 11980 1) are overlaid with our best fit broken power 
laws for each. Note the the mass-radius fit has two breaks as opposed to only one in the mass-luminosity fit. 
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Fig. 10. — The mass-temperature relation, again using data f rom lHarmaned 119881) and lPoppeJ 119801) . Overlaid onto the plot are our calculated relationship, 
as well as the analytical relation described bv lHarmaned ) 1 98^ . 
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Fig. 1 1. — Our calculated U-, B-, V-, and R-band magnitudes as a function of mass, compared to observed stellar magnitudes in those bands ICoxll2000h . 
Below each is a plot of the residuals in that band. 
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Fig. 12. — Our calculated I-, J-, and K-band magnitudes as of function of mass, compared to observed stellar magnitudes in those bands )Coxll2000l) . Below 
each is a plot of the residuals in that band. 
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Fig . 13. — A compaiison of our analytic formulation of the window probability versus using an exact formulation to calculate the average window probability 
as a function of total observing nights averaged over three different orbital period ranges assuming a linear distribution of periods. The left side uses a two transit 
minimum, while the right requires three transits. The top two plots - a) and b) - use the idealized conditions of 8 hour nights and no observing losses due to 
weather The middle plots - c) and d) - also use 8 hour nights, but now 1/3 of the total number of nights are lost due to bad weather. Finally, the bottom two plots 
- e) and f) - show the effects of shorter 4 hour nights but no losses due to weather. 
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Fig. 1 4. — a comparison of the PDMFs proposed by lReid et alj 420021) and lChabrie3 <20O3f ). IReid et alj 120021) uses a power law (equation |C II ) below a solar 
mass, while lChabrieJ 420031) uses a log-normal function (equation IC2I ). 
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Fig. 15. — The effects of various scale lengths on the predicted number of detections in a 6° X 6° field. The left shows detections as a function of Galactic 
longitude along the Galactic plane, as well as the difference in the number of detections from the assumption of a scale height of 2.5 kpc. The right plots 
detections as a function of Galactic latitude for three different Galactic longitudes. Note the lack of effect of the scale length at / = 90° , as well as the convergence 
of all three at fo = 90° , as is to be expected. 
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Fig. 16. — The total number of detections in a 6° X 6° field as a function of Galactic latitude for a variety of scale height relations. Also shown are the 
number of detections aro und stars w ith My < 5 and 5 < My , as well as the difference in t he predicted number of det ections using the constant 290 pc and Allen 
Astrophysical Quantities jCoxl2000l) scale height relations relative to the predictions using lBahcall & Soneiral 119801) . 



